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Abstract 

The gaussian ensemble and its extended version theoretically play the important role of interpolating ensembles between the 
microcanonical and the canonical ensembles. Here, the thermodynamic properties yielded by the extended gaussian ensemble 
(EGE) for the Blume-Capel (BC) model with infinite-range interactions are analyzed. This model presents different predictions 
for the first-order phase transition line according to the microcanonical and canonical ensembles. From the EGE approach, 
we explicitly work out the analytical microcanonical solution. Moreover, the general EGE solution allows one to illustrate in 
details how the stable microcanonical states are continuously recovered as the gaussian parameter 7 is increased. We found 
out that it is not necessary to take the theoretically expected limit 7 — > 00 to recover the microcanonical states in the region 
between the canonical and microcanonical tricritical points of the phase diagram. By analyzing the entropy as a function of 
the magnetization we realize the existence of unaccessible magnetic states as the energy is lowered, leading to a breaking of 
ergodicity. 

PACS numbers: 05.20.Gg, 05.50.+q, 05.70.Fh, 65.40.Gd 
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I. INTRODUCTION 

The canonical and grand-canonical ensembles approx- 
imate the microcanonical ensemble in the limit of in- 
finitely large number of particles, where surface effects 
and fluctuations can be disregarded with respect to the 
bulk mean values [l], 0). However, if the system sizes are 
not large enough compared to the range of interactions, 
or even in the presence of long-range forces, this inher- 
ent expectation changes dramatically. Although such 
nonextensive systems can be appropriately described by 
models in a volume-dependent scaling manner Q, the 
non-additive character still remains. As a matter of fact, 
the lack of additivity can be noticed for astrophysical 
objects, where gravitational interaction is responsible 
for a nonnegligible contribution from particles at large 
distances Thus, one realizes that most of the systems 
in nature can be encompassed in a class that can be 
designated non-additive for what concerns energy and 
entropy. The existence of a fundamental ensemble, the 
microcanonical one, seems to meet a consensus, while 
others, in particular the canonical and grand-canonical 
ones, are taken as its approximations (USUI]- In cases 
where full consistency of statistical ensembles holds for 
systems that undergo phase transitions, it is found that 
finite size scaling relations still place the microcanonical 
approach as the fundamental one [7J. 

There are many examples of systems whose equilibrium 
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properties are not equivalent in both microcanonical and 
canonical ensembles. Differences in the thermodynamic 
features have been verified analytically for systems with 
long-range interactions fcMTTl]. These examples show 
that the nonequivalcnce appears where the canonical 
ensemble presents a phase diagram with a first-order 
transition line. Actually, necessary and sufficient con- 
ditions for equivalence of ensembles can be formally 
stated [12( 1 . Thus, apart from the expected difference 
in the intermediate values of extensive thermodynamic 
quantities when one works with finite systems, the non- 
additive property also sets striking differences in the 
thermodynamic limit, leading to different phase diagrams 
[l3L Il4j |. Such nonequivalence has its counterpart in 
the nonconcavity of the entropy as a function of energy, 
S = S(E) [1, [ill EH- This may result in uncommon 
features at first-order phase transitions like temperature 
discontinuity and negative specific heat in the micro- 
canonical ensemble (3 [lH, [l5|, ll7HT9l |. In turn, the 
thermodynamic temperature (3 — 1/T(E) = dS/dE, (we 
take the Boltzmann constant ks = 1) is not a monotonic 
function of the energy and the equilibrium value E{j3) 
may be a multivalued function of /3 [H, Is [H| . 

An alternative ensemble, the gaussian ensemble [20I — 
l24l | was introduced to deal with systems that exchange 
energy with a finite reservoir. This contrasts with the 
canonical ensemble, where the system is in thermal 
contact with a huge heat reservoir and the energy ex- 
change is controlled by the temperature of the reservoir, 
which defines the average energy of the system. On 
the other hand, in the limit of no energy exchange with 
the reservoir, the system is isolated and thus has fixed 
energy. This is the microcanonical point of view, whose 
experimental situation resembles a system in contact 
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with a fictitious reservoir of extremely small size where 
the energy exchange can be disregarded. The gaussian 
ensemble has also been described as a regularization pro- 
cedure for the microcanonical ensemble [25| . Later on, 
Johal et al. [26j redefine the assumptions characterizing 
the gaussian ensemble to describe the thermodynamic 
properties of a system also in contact with a finite 
reservoir. This led to an extended version of the former 
gaussian ensemble. The extended gaussian ensemble 
(EGE) presents a smooth interpolation between its lim- 
iting behaviors, corresponding to the microcanonical and 
canonical ensembles. 

This work explores the EGE as a working ensemble 
for a system where the ensemble nonequivalence has 
been demonstrated, the Blume-Capel model. We ex- 
plicitly work out the analytical microcanonical solution 
from this ensemble. By means of the general EGE 
solution we are able to illustrate in details how the 
stable microcanonical states are continuously recovered 
as the gaussian parameter 7 is increased. We investigate 
the EGE behavior in the region of the phase diagram 
where one observes a first-order phase transition line 
in the canonical ensemble but of second-order type in 
the microcanonical description. Then, we point out 
how EGE identifies the canonical and microcanonical 
tricritical points. Moreover, we call attention to the 
broken ergodicity found in this model. 

The EGE formulation encompasses a natural extension 
of Statistical Mechanics, to include non-additive systems. 
The relation between EGE and Tsallis statistics has been 
described in Ref. (2(1 with the Tsallis parameter 
q being related to the parameter 7 in the EGE. The 
theoretical background characterizing this ensemble is 
presented in Section 2, where we briefly review some 
thermodynamic relations that are 7 dependent. The 
EGE solution of the mean-field BC model is carried out in 
Sec. 3 and is confronted at the thermodynamic level with 
the usual solutions in the canonical and microcanonical 
ensembles. The main conclusions about the effectiveness 
of the EGE in determining thermodynamic properties are 
summarized in Sec. 4. 



II. EXTENDED GAUSSIAN ENSEMBLE 

The canonical ensemble describes thermal properties 
of a system in thermal equilibrium with a heat reservoir. 
A new insight has been obtained when the reservoir is 
finite and possibly small. To this end, let a be a system 
with energy E and entropy S, and b a reservoir with 
energy E b and entropy S b , which exchanges energy with 
a. As a consequence, the energy of the system is allowed 
to fluctuate. Both systems form an isolated system 
with total energy E t = E + E b and total entropy St- 
Equilibrium is reached when the total entropy S t (E) is 
a maximum. The system itself and its heat bath can 
be considered subsystems of an isolated system where 
E fluctuates around its mean value U. Thus, for fixed 



external parameters like total energy E t and number 
of particles in the system, the most probable energy U 
is such that the expansion of the reservoir entropy S& 
around its equilibrium value E t ~U can be written up to 
the second order as 

S b (E b ) = S b (E t -U)+ (^) (U - E) 
\dE b J Et _ v 



1 fd 2 S t 



2 V dE 2 j Et _ v 



(U -E) 2 + 



(1) 



Because the derivatives depend on the reservoir thermo- 
dynamic properties, one defines [26| 



and 



dS b 
dE b 



d 2 S b 
dE 2 
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E t -U 
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E t -U 



In the case of an infinite reservoir, one would be working 
with the canonical ensemble and a would thus be iden- 
tified with the inverse thermodynamic temperature, a = 
l/T. This is because in the canonical ensemble approach 
the temperature T of the reservoir is a fixed parameter 
that determines the mean energy of the system. The 
effect of an infinite reservoir with constant temperature 
yields 7 = and vanishing higher-order derivatives in 
Eq.([T]); otherwise, those terms should be taken into 
account. 

The EGE is defined by the condition 7 7^ and 
probability density 



PjA E ) = 



p(E)e 



-aE—y(E-Uy 



Z 1 (U, a) 



(4) 



where Z 1 {U,a) stands for the normalization constant, 
which is the corresponding partition function in EGE 
[26l [28l I29I ] , with density of states p(E) and parameters 
7, a, and the dependent one U — U(a, 7). Actually, 
the extended gaussian ensemble is a particular case in 
a class of general functions g(E) [28|, [29|; the quadratic 
form g(E) — -f(E — U) 2 is just a convenient choice. The 
probability density in Eq. Q can be used to write the 
average energy of the system, 



U = j EP~ ha {E)dE. 



(5) 



Let us also introduce the extended thermodynamic 
potential analogous to the one in the canonical approach, 
$ 7 ([/, a) = — lnZ 7 (Z7, a). From here, the derivative at 
fixed value 7 can be obtained, 



d<S>-, 
da 



U, 



(6) 
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which parallels that of the usual canonical approach. The 
average energy U can be found self-consistently by means 
of Eq. ([5|) , which recovers the usual canonical ensemble 
result for 7 = 0, or from Eq. ©, as describing the 
equilibrium average energy with fixed parameters 7 and 
a. In this paper we follow a kind of inverse problem, 
U will be set as an input parameter that must, in con- 
junction with the variational problem of minimization of 
the extended thermodynamic potential, satisfy stability 
conditions for some (unknown) temperature 1/a, which 
is U dependent. The extended heat capacity has also 
been introduced [H H , C 7 = -a 2 {dU/da) 1 . 

The usual canonical ensemble deals with homogeneous 
configurations in equilibrium as a function of intensive 
variables like temperature. The canonical averages 
always produce smooth distributions of mean values, 
as in the case of heat capacity, when averaged over 
fluctuations. In contrast to the canonical heat capacity, 
the extended heat capacity may present negative values 
when 7 > 0. Negative values of C 1 {U) require that 
(dU /<9a) 7 > 0. Thus, van der Waals loops can be seen in 
this formalism. This sort of behavior has been observed 
in typical caloric curves, temperature versus mean ener- 
gies, for systems with thermodynamic first-order phase 
transition, a forbidden phenomenon in the canonical 
picture. These features are illustrated in Fig. 2 for the 
Blume-Capel model. Thus, the standard homogeneous 
thermodynamics given by the canonical ensemble is not 
suited to describe first-order phase transitions. On the 
other hand, the stability condition of a system is related 
to the homogeneous temperature that defines thermal 
equilibrium with the huge reservoir. The EGE includes 
the possibility of a rather small heat bath, thus allowing 
for the appearance of inhomogeneous configurations in 
the system for finite 7, which results into a weakened 
version for the constraint of constant energy that defines 
the microcanonical ensemble. 

The extended entropy can be obtained by the 
Legendre-Fenchel (LF) transform of the extended canon- 
ical thermodynamic potential $ 7 (a) as (2(| H^, [2l| 

From this transform and Eq. it follows that a(U) — 
dS 7 /dU. 

Notice that the above relations recover the canonical 
results in the limit 7—7-0. In this case, one has 
the standard Legendre transform S(U) = (3U — $(/3), 
where $(/3) = KrOy-^o^-yCo! = (3) corresponds to the 
canonical potential and U is the equilibrium mean energy, 
U = d&(/3)/dj3. It is well known that the standard 
Legendre transform of $(/?) always produces a concave 
function of U. Therefore, nonequi valence between mi- 
crocanonical and canonical ensembles appears when the 
microcanonical entropy is a nonconcave function of U in 
some energy range, as shown in Fig. 2b. In that case of 
nonequivalence, S(U) can be named as just the canonical 
entropy: S can = /3U clin (J3) + In Z can (/3). On the other 



hand, the limit 7 — > 00 corresponds to th e microcanonical 
case. This can be seen as lim 7 _ i . 00 y / 7r/7Z 7 ([/, a) through 
the use of the Dirac's delta sequence in the gaussian form 
[25j]. For finite 7 one obtains an intermediary thermal 
description between the known limiting ensembles. 



III. EXTENDED GAUSSIAN SOLUTION OF 
THE MEAN-FIELD BC MODEL 

The Blume-Capel model is a spin-1 Ising model |30l 
l3~i| and was introduced to describe phase separation 
in magnetic systems. It is a particular case of the 
Blume-Emery-Griffiths model [32| aimed at describing 
the critical behavior of He 3 -He 4 mixtures with different 
concentrations. Here we consider its mean field version, 

N / N \ 2 

H(S) = Aj2sf-^ K> , (8) 

i=l \i=l / 

where Si = {0,±1}. The couplings J > and A are 
the exchange and crystal-field interactions, respectively. 
The BC model represents a simple generalization of the 
spin-1/2 Ising model, but with a rich phase diagram in 
the (A/ J, T/J) plane. It exhibits a first-order transition 
line, tricritical point, and a second-order transition line. 

The critical properties of the BC model can be de- 
termined analytically in both the microcanonical Q 
and canonical ensembles [H, [33j]. It has been demon- 
strated that these ensembles do not yield the same 
phase diagram for the first-order critical line Q. The 
canonical tricritical point occurs at (A/ J, T/J) = (~ 
0.46209812,1/3), which gives origin to the first-order 
transition line for larger values of A/ J. The mi- 
crocanonical solution identifies the tricritical point at 
(A/ J, T/J) ~ (0.46240788, 0.33034383). 

Here, it is useful to introduce the order parameters 
magnetization M = Xh=i &i = — N- and its second 
moment, the quadrupole moment Q = YliLi &i = + 
N-, where N + and _/V„ are, respectively, the number of 
sites with up and down spins. If No is defined as the 
total number of zero spins, then N = N + + iV_ + A*o is 
the total number of spins in the system. 

The extended gaussian partition function 

Z,(U, a) = Y, e-^W-iMW-vp f (Q) 

{S} 

can be analytically solved in terms of its order param- 
eters M and Q. To this end, the so-called Hubbard- 
Stratonovich (HS) transformation 

+00 




d ye - y2/2b - lxv , (10) 



is applied to the gaussian term in Eq. © with the choices 
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b = 2j and x — H(S) — U. It turns out that 



{S} 



(11) 



+ 00 



dye -^-(w+^(AT,t 1 Sf-U) e (iv+ a )(y/ i L r j:^ 1 S l ) 2 _ 



(12) 

By making use of another HS transformation, the 
extended gaussian-partition function becomes 



^^tt^E 



-all 



{S} 



+00 

x dy 



1/2 



*V + a \ ' e -3^-(iy+a)(AEjLiS?-E0 



+ OO 

J dze^ {tv+a)z2+2{tv+a)z ^ r S^=^ Si . (13) 



Since 5j = {0, ±1}, it follows that 

^ e -(iy+°<) A £i S*+2(iy+a)z Si 



{S} 



= 1 + e 



-(iy+ 



a)A ^ e 2(iy 



+a)z 



-2{iy+a)z -> 



N 



(14) 



AT N-N 

E E 

AT =Q AT + =0 



AH 



N \N + \N_\ 



,-(iy+a)A(N-N ) 



-2(iy+a)zJ^{N-N a -N + ) 2{iy+a)zJ^N + 



xg 1 -I- v 2H"-"«-"+/ e *r»T"i' v 51"+ ^ (15) 

where the last result is obtained by applying the binomial 
expansion twice to the result in Eq. (|T4"]) . Now, placing 
the order parameters M and Q in Eq. (|15p and inserting 
this result into Eq. (TIB"]) , one obtains 

1 N N ~ No yvi 

V / 7V o =0iV + =0 U + 



+00 



1/2 



x / dy [ l l±^\ e -*H<*+«>«-*.' / 



(16) 



This expression can be integrated by gaussian formulas 
to produce 



N N-N 



Z 1 (U,a)= E 



AH 



x e 



iV !./V+!iV_! 

-Q(AQ-iM 2 )- 7 (AQ-t/-iM 2 ) 2 



(17) 



The solution for this ensemble notoriously brings forth 
the counting factor for the number of microscopic states 
corresponding to the macrostate defined by M and Q. 
These order parameters indeed define the energy E of a 
configuration given by the Hamiltonian in Eq. ([8jl, 



E = AQ-±M*. 



(18) 



Hence, it is convenient to write explicitly the extended 
partition function as a function of those order parame- 
ters, 



Zj(U, a) 



N M=Q 

E E 

Q=0 M=-C 



AH 



(N-Q)l [| (Q + M)]![|(Q-M)]! 



x e - Q (AQ-^)- 7 (AQ-^-C/) 2 



(19) 



Before studying the thermodynamic features presented 
by this ensemble as a function of finite 7, it is important 
to show explicitly the limiting microcanonical behavior 
of this ensemble. 



A. Microcanonical limit and negative response 
functions 



Let us firstly explore the limit 7 —¥ 00 to obtain the 
microcanonical ensemble. Since it is required that the 
extended partition function is well-behaved in this limit, 
the sum in Q must converge to a dominant value for some 
Q such that AQ - JM 2 /2N - U = 0. This is nothing 
else than the microcanonical constraint on the energy E. 

Next, the thermodynamic limit A^ — > 00 is studied. 
Here, it is convenient to work with the intensive quan- 
tities q — Q/N and m — M/N. Let us also define 
K = J/2 A and e = U/AN as in Ref. [8]. Equation 
(fl~8]| now reads e = q — Km 2 and becomes a constraint 
equation for the average energy e as 7 — » 00. 

For large N, one can evaluate the microcanonical parti- 
tion function Z{e,a), where Z(e,a) — lim 7 _ ) . 00 Z 7 (e,a), 
as a variational problem. To this end, we consider the 
saddle point solution, Z(e,a) ~ e —Nip(e,a,m) ^ wnere m \ s 
such that the thermodynamic potential 




+ ln(l-g) 



(20) 

is minimized for each average energy per site e. To obtain 
tp(e, a, to), we also applied the Stirling approximation for 
large N to Z(e,a) in Eq. (TT9"|) . The above expression, 
Eq. (f2"U|). was kept as a function of q and to, to 
recognize that the term inside the brackets is the correct 
microcanonical entropy — s m j cro (e, to) obtained in [§[ as 
a function of the parameter mean energy e and mean 
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FIG. 1: Entropy s m i cro (e, m) for some values of e with 
A/J — 0.462407. Gaps in the magnetization correspond to 
unaccessible states. 
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magnetization m. The entropy is an even function of m 
and a nonconcave function of the independent variables 
m and e, as respectively shown in Fig. 1 and 2b. This 
fact has striking consequences for the response functions 
specific heat and specific susceptibility. 

Here, we remark that the microcanonical entropy is not 
always an analytic function. As a consequence, gaps may 
develop in the magnetization for some values of s as illus- 
trated in Figure 1 for the coupling A/J = 0.462407. This 
means that the system presents ranges of disconnected 
magnetization as function of e and A/J. Thus, we cannot 
move continuously from one domain of magnetization 
to any other, leading to the so-called microcanonical 
ergodicity breaking (9j , which is not related to any phase 
transition. The condition for unaccessible magnetization 
states can be easily determined from the expression for 
the entropy. Thus, one finds that those gaps start at 
e = A/2 J and increase as the energy e is lowered. 

Now, back to the EGE approach, it is worthy of 
mention that we are not evaluating a Laplace integral 
of the usual canonical partition function: the extended 
thermodynamic potential per site </? 7 results from a 
modified partition function 

<^ 7 (£,o;) = — lim — lnZ 7 (e,a), (21) 

which transforms nonequilibrium states of the canonical 
ensemble into equilibrium states of the extended ensem- 
ble. Here, the dependence of <p-y on m has been omitted 
because we are already assuming that the minimization 
in m has been accomplished. As emphasized, the 
nonequivalence of ensembles (microcanonical and canon- 
ical) is a consequence of the anomalous behavior of the 
microcanonical entropy characterized by the existence 
of convex parts in s m i cro (e,m). The nonconcavity of 
the entropy function means that the system contains 
several energy-dependent equilibrium states, revealed in 
the microcanonical ensemble, which do not have their 



FIG. 2: Microcanonical behavior of the BC model with 
A/J = 0.462407. (a) Microcanonical temperature as a 
function of the average energy e. The horizontal dashed 
line corresponds to the canonical critical temperature of the 
transition, (b) The shifted microcanonical entropy s(e) = 
s m icro(e) — (A+Be). The subtraction is performed to visualize 
the nonconcavity of the entropy in relation to the linear 
function joining s m i cro (e a ) to s m icro(£i>). (c) Specific heat c(e). 
It presents two poles located by the zeros of the determinant 
d(e,m), where m stands for the values that maximize the 
entropy at e. Those poles can also be read from T(e) behavior 
in (a). c(e) becomes negative in between those poles, (d) 
Specific susceptibility xi e )- It presents two poles, again 
placed at the zeros of d(e, m) and becomes negative between 
them, (e) Behavior of the determinant d(e, m) as a function 
of e. The vertical dashed lines signal the zeros of d(e,m). 



counterpart in the temperature-dependent equilibrium 
states in the canonical description. Thus, the new term 
in 7 turns such points s(T) into equilibrium points in the 
extended ensemble [5^, [29|, [U, [35[ . The usual thermo- 
dynamic potential ip(e, a) is given by the minimization 
procedure 

ip(e, a) = mini^fe, a, m) , (22) 

m 

where the dependence of ip on e is always kept to show 
that <p(a) is calculated at the equilibrium value that 
minimizes this potential. 

In the present ensemble, the LF transform Q of 
<y9 7 (e, a, m), where e, a and m are independent variables, 
produces the correct s m i cro (e) as follows, 

s m icro(£) = min max { lim s 7 (e,a,m)}, (23) 

a m 7— >oo 

where lim 7 _>. 00 s 7 stands for eaA — tp(e,a,m) in this 
model. From this result one recovers the known ther- 
modynamic behavior. Figure 2 contains our calculations 
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for the microcanonical temperature T(e), shifted entropy 
s(e) = s m i cro (e) — (A + Be), specific heat, susceptibility 
and determinant of the curvature of s m i cr o, as a function 
of e for A/J = 0.462407. This value of A/J is in the 
canonical first-order phase transition region but it is in 
the microcanonical second-order phase transition region. 
Since e = U /AN, one obtains 



1 



1 ds T 



T(e) 



de 



(24) 



The horizontal dashed line in Fig. 2a indicates the 
temperature T can ~ 0.330666 obtained by canonical 
methods [H, [33|. It connects the point T(e a ) to T(eb), 
where e a ~ 0.328959 and ~ 0.330646, are read from 
Fig. 2a. The width 5e = 0.001687 is the specific 
latent heat of the first-order phase transition seen in the 
canonical ensemble. 

In Fig. 2b it is shown s m i C ro(e) shifted by the canonical 
entropy s(e) = A + Be, where A ~ 0.401447 and B ~ 
1.398397 are such that s(e a ) — s m i cro (£ a ) an d s(eb) — 
Smicro (e&)- This subtraction allows one to highlight the 
so called convex intruder in the specific entropy @, @ ■ 

The point c, as signaled in Fig. 2a, corresponds to 
the energy e c where occurs the minimum of the shifted 
entropy Points d and e in Fig. 2b signal the energy 
range (ed,e e ) where the entropy is nonconcave, ea — 
0.3297040 and e e ~ 0.3303532. In Fig. 2a, we have the 
corresponding temperature T(ed) — 0.33074967 as the 
maximum temperature in that energy range. 

Figure 2c shows the specific heat 



c(e) 



de 
dT(e) 



T 2 d{e,m) 



(25) 



where d(e,m) is the determinant of the curvature of 

Smicro (^0 ? 



die, m) — det 



Sine $mm 



(26) 



where the notations s 



and 



refer respectively 



to the second derivatives d 2 s m i cro /de 2 , d 2 s m i CIO /dedm 
and d 2 s m i cro /dm 2 . This determinant addresses the 
stability conditions around the stationary points m and 

e 

Figure 2d presents the corresponding magnetic suscep- 
tibility 



3 (e, m) 



r/A 2 



d(e, m) 



(27) 



The nonconcavity of the microcanonical entropy in e 
and m renders a negative region for the specific heat 
and magnetic susceptibility. Finally, Fig. 2e depicts 
the behavior of the determinant d{e, to) as a function 
of e with to evaluated at the microcanonical equilibrium 
condition. The zeros of this determinant indicate the 
region where the response functions attain negative val- 
ues. Their locations are represented by vertical dashed 



lines. Those negative values for the canonical observables 
happens inside the convex region related to the phase 
separation in the first-order thermodynamic transition. 



B. 
tial 



Finite 7 and extended thermodynamic poten- 



For finite 7, one obtains different equilibrium prop- 
erties. As we are going to show, the full equivalence 
with the ones in the microcanonical ensemble is achieved 
for finite 7 only for A/J between the canonical and the 
microcanonical tricritical points. On the microcanonical 
first-order transition line, one needs 7 — > 00 for such full 
recovery of the microcanonical results. 

The analytical solution for the extended thermody- 
namic potential is analogously obtained following the 
procedure leading to Eq. (|2TJ|) . 

p 7 (e,ar,m,g) = 



gin 



n l \ m , / q + m , 



+7A 2 (g - Km 2 - ef + aA(q - Km 2 



(28) 



The basic problem concerning nonequivalent ensembles 
is that the true s m i cro cannot be obtained as an LF 
transform of the free energy (p(j3). Here, the application 
of the extended LF transform to tpy yields the extended 
entropy s 7 , which can be read from Eq. (|28|) . s 7 = 
aA(q — Km 2 ) — ip^. This entropy is now a concave 
function of e. The extended inverse temperature reads 
a = (dsj/de)/ A. It characterizes stationary points 
analogous to the physical inverse temperature and is 7 
dependent. 

Now, let us evaluate the equilibrium points of 
f-y(e,a,m,q). Notice that the microcanonical constraint 
for the specific quantities is not enforced here, the 
variables e, to and q are treated as independent variables. 
The linear term in 7 can be seen as a constrained 
equation, leading to the microcanonical ensemble only 
for 7 — > 00. 

It was verified that all solutions of dip^/da = eA, as 
expressed by Eq. and dipj/dm = 0, dip 1 /dq = 0, 
for a fixed e, are only the ones given, for example, for 
T(e) in Fig. 2a. However, those solutions are not stable 
for all 7. Since the analytical expression cp~ comes from 
the saddle-point approximation in Eq. (|19[) , one needs to 
study the stability of those EGE solutions as a function 
of to and q. To this end, the determinant of the Hessian 
matrix, 



d(m, q) = det 



dm 2 
"Vr 



dmdq 



dqdm Bq 2 



is analyzed in the T versus e plane as a function of 7. This 
amounts to exploring which points {m, q} minimize tpj 
for fixed T and e, and satisfy the condition d(m, q) > 0. 
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condition on 7, 



FIG. 3: EGE temperatures for some values of 7 with A/J 
0.462407. 



y 




0.328 0.3285 



0.331 



FIG. 4: 7 = — l/2AT 2 (e)c(e) presents positive values for e € 
(Ed,£e): energy range of nonconcave entropy. Negative values 
of 7 occur at energies out of this range. Here A/J = 0.462407. 



Figure 3 shows the lines of stable points for different 
values of 7. Notice that the figure for 7 = corresponds 
to the canonical results, but we have not included the 
Maxwell construction. All presented states are the stable 
ones in the canonical ensemble. This procedure selects 
solutions T(e) for energies where the entropy is a concave 
function. The gap in e corresponds to the region where 
one observes negative values in the specific heat and 
in the susceptibility. As 7 increases, one recovers the 
microcanonical solution. In fact, for sufficiently large 
7, s 7 (e) becomes entirely concave and continuous on 
ee(e d ,e e ) [H,|1I, 



ds 



< 0. 



(29) 



The addition of the term in 7 to the usual Legendre 
transform changes the energy range where the noncon- 
cavity of the canonical entropy is observed. How this 
energy range is reduced as 7 increases can be easily 
evaluated from Eq. (|2"9"|) . This implies the following 



7 > 



2AT 2 (e)c(e) ' 



(30) 



But in view of the specific heat c(e) < for s € 
(£(i,£ e ), one obtains 7 > in this range, as exhibited 
in the inset of Fig. 4. Figure 4 shows the behavior 
of 7 = — l/2AT 2 (e)c(e) for energies out of that range, 
too. This figure highlights the minimum value of 7 to 
achieve equivalence with the microcanonical ensemble 
for s £ (ed,s e ). The full equivalence in this energy 
range is reached when 7 ~ 4950 for the example with 
coupling A/J = 0.462407. Negative values for 7 have 
been considered in [27j to enhance Monte Carlo sampling. 
Here, a negative 7 converts microcanonical stable states 
at e to unstable ones when e < e a or e > £{,. Figure 5 
shows, for all values of A/J between the canonical and 
microcanonical tricritical points, the minimum 7 needed 
to recover the exact microcanonical solution. From 
the canonical approach, a first-order phase transition 
starts at A/J ~ 0.46209812, but from a microcanonical 
analysis the true first-order transition starts at A/J ~ 
0.46240788. The EGE approach distinguishes those 
transition regions presenting finite values for 7, to recover 
the full thermodynamic features of this model when A/J 
is between those tricritical points. 



10000,= 




0.01 



0.4621 



0.4622 



0.4623 



0.4624 



A/J 



FIG. 5: Minimum 7 needed to recover the full microcanonical 
solution. The vertical dashed lines signal the canonical 
(A/J ~ 0.46209812) and microcanonical (A/J ~ 0.46240788) 
tricritical points. 



IV. CONCLUSIONS 

In conclusion, the analysis of the BC model shows how 
the stable states present in the microcanonical approach 
but not found in the canonical one can be obtained from 
EGE. This approach leads to analytical expressions for 
the extended free energy and entropy in a simple way, 
and quantifies the nonequivalence of ensembles between 
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the tricritical points. The EGE formulation exhibits 
negative specific heat, like the microcanonical one, in the 
canonical first-order phase transition region. This also 
happens between the tricritical points where 7 is finite. 
As a consequence, this remark may open a way of finding 
tricritical points in systems where analytical solutions 
can not be obtained. Thus, an appropriated Monte Carlo 
method based on EGE should be preferable than the 
standard one where sampling relies on the Boltzmann 
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